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ABSTRACT 

We describe the production and verification of sky maps of the five SPIRE fields observed as 
part of the Herschel Multi-tiered Extragalactic Survey (HerMES) during the Science Demon- 
stration Phase (SDP) of the Herschel mission. We have implemented an iterative map-making 
algorithm (SHIM; The SPIRE-HerMES Iterative Mapper) to produce high fidelity maps that 

^-^ preserve extended diffuse emission on the sky while exploiting the repeated observations of 

the same region of the sky with many detectors in multiple scan directions to minimize resid- 
ual instrument noise. We specify here the SHIM algorithm and outline the various tests that 

^-H were performed to determine and characterize the quality of the maps and verify that the 

!• astrometry, point source flux and power on all relevant angular scales meets the needs of 

. the HerMES science goals. These include multiple jackknife tests, determination of the map 

transfer function and detailed examination of the power spectra of both sky and jackknife 
^ maps. The map transfer function is approximately unity on scales from one arcminute to one 

degree. Final maps (vl.O), including multiple jackknives, as well as the SHIM pipeline, have 
been used by the HerMES team for the production of SDP papers. 
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by European-led Principal Investigator consortia and with important partic- 
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1 INTRODUCTION 

The Spectral and Photometric Imaging Receiver (SPIRE; Grif- 
|fin et al.||2010|) o n board the Herschel Space Observatory (Pil- 
bratt et al.||2010) has opened a new window on the Universe at 
far-infrared/submillimetre (FIR/submm) wavelengths. The spec- 
tacular sensitivity of SPIRE detectors combined with HerscheVs 
3.5 m aperture, currently the largest telescope in space, allows as- 
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tronomers to observe the FIR/submm Universe with unprecedented 
efficiency. During the Science Demonstration Phase (SDP) of the 
Herschel mission, five fields of various size and depth were ob- 
served as part of the Herschel Multi-tiered Extragalactic Survey 
(HeMES, |Qlrver et al.|2010fl . 

The SPIRE photometer focal plane arrays (FPAs) consist of 
139, 88 and 43 feedhorn-coupled spider-web bolometers which 
measure temperature changes resulting from the absorption of in- 
cident radiation in three bands centered on 250, 350 and 500 |im 
respect ively. The Herschel Interactive Processing Environment 
(hip^J Ott 20 10) data processing pipeline, which converts raw 
bolometer data to physical units and removes artefacts like cos- 
mic rays and large-scale signal variations due to the temperature 
changes in the FPAs are removed is described in detail in Dowell et 
|aLl (2010a> . 

This pipeline is broken into several stages, including the con- 
version of raw telescope data into engineering units (i.e. Volts), the 
conversion of voltage into flux density units, corrections for elec- 
trical, optical, thermal and astrophysical artefacts and, finally, map- 
ping of the time-ordered data into sky images. Details of the hipe 
settings used for the current processing are given in Section |3.1| 
hipe works on detector timelines which consist of the detector sam- 
ples for a single scan, at constant scan speed, across the map area. 
Currently, the default map-maker in hipe removes a median from 
each detector timeline and then simply bins the measured fluxes 
into map pixels. For the majority of the HerMES Science Demon- 
stration Phase (SDP) science goals, particularly those studies that 
focus on the detection and characterization of individual sources, 
the hipe 'naive' maps meet all the requirements. However, science 
goals that depend on the statistical power in the map to measure 
large-scale structure and constrain source count models, are partic- 
ularly sensitive to artefacts and/or residual power (or lack thereof) 
in sky maps due to the processing. 

The SPIRE FPAs are held at « 0.3K and the temperature fluc- 
tuations due to sources on the sky are generally much smaller than 
slow temperature drifts in the focal plane structure itself. For this 
reason, a single, raw scan of the sky is dominated by low frequency 
noise, correlated among detectors, and commonly referred to as 1// 
noise. The SPIRE instrument tracks these drifts using both thermis- 
tors on the focal plane and dark bolometers that are not exposed to 
the sky. These additional data allow removal of the \/f noise by 
the hipe pipeline. The drifts are removed so completely by the stan- 
dard pipeline that, as stated above, for studies of individual sources, 
there is no appreciable gain in achieved sensitivity in going from a 
simple binning mapper to a more sophisticated method. For sta- 
tistical studies, however, any subtle phase difference in thermis- 
tor/bolometer response, an uncoupled bolometer drift, or an error 
in interpreting the thermistor coupling to the bolometer signal can 
introduce or leave false signal in the maps that, while difficult to de- 
tect, can affect science results. Additionally, the process of binning 
time-ordered data into sky maps must be optimized to accurately 
reproduce optical power on all relevant angular scales and must be 
characterized for potential map making artefacts that can bias sta- 
tistical interpretations of map data. 

This paper describes the production and characterization of a 
set of high-fidelity maps, now in use by the HerMES team, pro- 
duced with a mapping algorithm designed to mitigate these poten- 
tial artefacts, should they exist in the post-HIPE data. 

In selecting a mapping algorithm for the HerMES SDP 



http://herschel.esac.esa.int/HiPE_download.shtml 



maps, it was anticipated, based on ground-based and balloon- 
borne FIR/submm/microwave experiments, that maximum like- 
lihood map-makers, such as MADmap (Cantalu po et al.| [2010) 
or SANEPIC jPatanchon eT1uH|2009] l would be needed. While 
MADmap has been implemented in hipe, it is not optimal for SPIRE 
data due to the large number of detectors with correlated noise. 
SANEPIC was designed for the data from BLAST, (the Balloon 
Borne Large Area Submillimeter Telescope, Patanc hon et al.|2 009) 
which flew a prototype of the SPIRE instrument on a stratospheric 
balloon and is an optimal map-maker made to specifically handle 
correlated noise among many detectors. Ultimately a SANEPIC- 
like optimal solution will be used to map the SPIRE data. How- 
ever, the environment at the second Lagrange point (L2) of the 
Sun-Earth system, 1.5 million km from Earth, is far more stable 
than any environment here on Earth, as well as having no atmo- 
sphere. Thus, it was determined early in the Herschel mission that, 
following temperature drift removal, the SPIRE data is free enough 
from \/f noise that most of the early SDP science goals could be 
achieved using a simple binning map-maker. For the more sensitive 
science goals, we have developed an interim map-maker based on 
the algorithm of |Fixsen et aTlpOOO] ). SHIM, the SPIRE-HerMES 
Iterative Mapper, exploits the redundant observations of each point 
on the sky with multiple detectors in cross-linked scan directions 
to determine both the sky signal and the weighting of the data sam- 
ples for binning into map pixels while simultaneously producing 
an accurate error map. Individual detector gains can also be de- 
termined, though in this instance the detector gains remain fixed. 
SHIM works with SPIRE-only and SPIRE-PACS^ Parallel modes. 
The full algorithm is described in Section [T2| 

Following implementation, verification that the maps are free 
of artefacts from both the low level pipeline processing of time or- 
dered data and of the mapping algorithm itself is required. In or- 
der to characterize the effects of the timeline processing and the 
mapper itself, a detailed examination of the astrometry, noise prop- 
erties, transfer function and other potential artefacts in the SHIM 
vl.O maps is described in Section [4] Future versions of the SHIM 
pipeline are expected to build upon what is described here. SHIM 
will be released as part of the HerMES initial data release (DR1), 
scheduled for the first quarter of 2011, via the HerMES website t 
and the algorithm will ultimately be coded into java by the NASA 
Herschel Science Center (NHSC) for inclusion in future versions 
of HIPE as a complement to the current naive mapper. 



2 DATA SETS 

The HerMES SDP fields are GOODS-N, Lockman-North and 
Lockman-SWIRE, FLS and the galaxy cluster Abell 2218. 
GOODS-N, Lockman-North and Abell 2218 were all observed in 
SPIRE Large Map mode at the nominal scan speed of 30 arcsec/s. 
The FLS field was observed in SPIRE-PACS Parallel mode, with 
a scan rate of 20 arcsec/s and the Lockman-SWIRE field was ob- 
served in SPIRE Large Map mode at the fast scan speed of 60 arc- 
sec/s. Full details of the observed field sizes, scan strategies and 
depths are presented in |Qliver et al.|p010a) . Sections |3^2| and |4~4| 
indicate that while the map depth, i.e. the number of detector sam- 
ples per map pixel, is affected by the scan speed, the mapping algo- 
rithm and the expected/detected map systematics are independent 



2 PACS: Photodetector Array Camera and Spectrometer, Poglitsc h" et al.| 
|2010| 
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of the observing mode among these fields. The SHIM vl.O maps of 
all five fields have been used by HerMES collaboration for several 
projects which rely on the statistical properties of the maps includ- 
ing (NguyerTet^^ 
and |Amblard et al. 1(2015) 7" 



3 MAP MAKING PIPELINE 

SHIM takes HiPE-processed timelines, calibrated in physical units, 
and injects them into an iterative mapper. Both phases of this pro- 
cessing are described below. 



3.1 Timeline Processing 

Raw telescope data were processed into calibrated timelines us- 
ing hipe version 3.0.1 108 with the latest (at the time of processing) 
calibration files and modules for removal of cosmic rays. This pro- 
cessing also includes a corrections, provided by the SPIRE Instru- 
ment Control Center, of a 79 ms time offset between the spacecraft 
and instrument clocks and a 2 Hz correction to the approximately 
300 kHz clock frequencies of the spacecraft and instrument clock ~| 
to fix a pointing offset seen in early observations. The details of the 
timeline processing are available in Dowell et al. (2010a) and in the 
SPIRE Observers Manual PSC-DOC-0789,|2010| ). In brief, once 
the raw data have been converted into detector voltages, several cor- 
rections to the timelines are made and the telescope pointing prod- 
uct is associated with each detector sample. The data are processed 
in the following order: first, detection and masking of cosmic ray 
glitches that affect all detectors ('frame hits') and then glitches that 
affect only individual detectors ('web hits'), in this case using the 
'simple' or 'sigma-kappa' deglitcher which considers timeline data 
points that are k x <t outliers to be glitches, where cr is the stan- 
dard deviation and k is chosen by the user. In this analysis we use 
k = 5. Data affected by cosmic rays are flagged and are always 
masked at the map-making stage. The signal timeline is interpo- 
lated over the affected samples to maintain a continuous timeline, 
but the interpolated data points are not included in the maps. The 
pipeline then corrects for the electrical filter response, converts sig- 
nal to flux density and removes temperature drifts based on the FPA 
thermistors. Finally, corrections are made for the bolometer time 
response before the signal timeline is merged with the telescope 
pointing product to produce sky coordinates for each detector at 
each sample. While modules exist for the correction of electrical 
and optical crosstalk, but currently make no correction. Calibrated 
timelines can then be fed into SHIM as described below. 



3.2 Mapping Algorithm 



Based on the algorithm described in |Fixsen et al.|p000| ), SHIM 
uses an iterative baseline removal and calculation of detector 
weights to produce 250, 350 and 500 \im maps with pixel sizes (6.0, 
8.3 and 12.0 arcsec) chosen to be approximately 1 /3 the FWHM of 
the beam in the relevant band. The data are dealt with in scans, i.e. 
one pass across the region of sky to be mapped. No timeline filter- 
ing is applied before the data are passed to the map-maker, except 
that the median of each scan is subtracted, which increases the effi- 
ciency with which the map-maker converges. Given sky brightness, 
M(x, y), the signal, S dsj , for detector, d, in scan, s, at time sample, 

3 These correction are now part in the standard hipe pipeline. 
© 2002 RAS, MNRAS 000,[l]|9] 



j, is modeled as a sum of the sky, an offset and measurement noise, 
i.e.: 



Sdsj = gd M(x dsj ,y dsj ) + p ds + N d5j , 



(1) 



where g d is the detector gain, and p ds is an order-rc polynomial 
baseline offset for detector, d, and scan, s. N dsj is the measurement 
noise, and (x ds j,y dsj ) encodes the position of each detector on the 
sky. The initial sky map M°(x ds j,y dsj ) is a 'naive' map in which the 
detector samples are simply binned into map pixels and the gains, 
g d are initialized to unity. If the detector gains, g d , are allowed to 
vary, the parameters g d and p n ds are each fit iteratively, with each 
iteration, i, consisting of two alternating steps. In the first step in 
each iteration, the g d are set to g l ~ l and are held fixed while the 
n + 1 parameters, p n sd , for each detector, d, and scan, s, are varied si- 
multaneously. Using a linear least squares minimization procedure 
the p\ are determined that minimize the variance of the timeline 
residual, R l dsj , based on the previous map, M l ~ l (x, y), 



Ksj = Sdsj ~ [§d M l -\x dsj ,y dsj ) + p2] • 



(2) 



In the second step, for each detector, d, g d is varied to minimize 
Rdsj, while the p n sd are fixed to p n sd . 

A new estimate of the sky is then constructed by taking the 
weighted mean of all samples that fall into each pixel: 



X w d S ( s dsj-p n d i)/g d 



M\x,y) ■■ 



dsj£(x,y) 



dsj e (x,y) 



(3) 



where the weight w l ds is the inverse variance of the timeline resid- 
ual, 



(4) 



and N is the number of samples in scan, s. Note that, while the 
weights are not free parameters, they are determined from the resid- 
uals as part of the iterative process. A noise map is created by prop- 
agating detector noise as estimated by the variance of the residuals: 



o-\x,y) = 



dsje(x,y) 



-1/2 



(5) 



The order n of the polynomial baseline is chosen based on the phys- 
ical size of the map where we use n = 1 for our smallest (9' x 9') 
cluster fields, n = 2 for intermediate (30' x 30') deep fields and 
scale up to n = 3 for the widest (2° x 2°) fields. To improve sta- 
bility of the algorithm, the iteration on which to begin fitting each 
component (or whether to fit it at all) can be specified. This in- 
cludes the choice of whether to hold the w ds fixed. The pipeline 
is typically run for several (5-10) iterations with the w ds fixed to 
unity and allowing w l ds to be determined by Eq. 4 in subsequent it- 
erations. In the maps presented here, the gains g d are fixed to unity. 
While the algorithm allows us to fit for the individual gains directly 
from the data this increase in complexity will be reserved for future 
versions of the maps. The detector gains for the SPIRE bolome- 
ters have been derived using dedicated observations of calibration 
sources and the g d should not vary significantly from unity at this 
stage. As a check, a single mapping run of SHIM on the GOODS-N 
field with g d allowed to vary produced g/s that differed from unity 
by less than a few percent. While this may lead to a marginal im- 
provement in the maps and will likely be implemented for future 
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Figure 1. A 0.7 deg x 1.1 deg sample of the FLS field displaying both the 
confused SPIRE point source background and a bright strand of Galactic 
cirrus. 



versions, the verification of this change is beyond the scope of the 
SDP phase. 

The method described in Eq. 5 of estimating the noise map 
differs from the hipe map-making pipeline, which uses the stan- 
dard deviation of all samples that fall into each pixel. Comparing 
the median noise in the error maps in both cases in the GOODS-N 
field, we find that they are not significantly different, indicating that 
the standard deviation provides a reasonable estimate of the over- 
all noise level in the map. The direct propagation of detector noise 
used here should give a better measure of noise variations across 
the map and is more robust when the number of samples per pixel 
is low and the standard deviation becomes a poor estimate of the 
statistical error in the measurement. While, in general, only a small 
fraction of pixels (< 0.1%) have less than 5 samples, in fast scan 
mode at 250 \im, where 7 dead or slow bolometers are masked, re- 
gions of the map that have had full scans masked due to temperature 
drift artefacts have pixels with as low as 2-3 samples. Splitting the 
data in half for the jackknife tests described in the following sec- 
tions then leads to occasional pixels with only a single sample, for 
which the standard deviation can not provide an error estimate. It 
should be noted that the FLS field, observed in SPIRE-PACS Paral- 
lel mode scanning at 20 arcsec/s in only a single repeat, has fewer 
low sample pixels than the two repeat observation of Lockman- 
SWIRE field observed in SPIRE Large Map mode scanning at 60 
arcsec/s just as one would expect given three times more samples 
per pixel scanning at 20 vs. 60 arcsec/s. 

High quality SHIM maps have been produced for all HerMES 
SDP fields. A 0.7 deg x 1.1 deg portion of the 250 [im FLS field is 
shown in Fig.[T] This particular field is chosen to deliberately illus- 
trate the range of structures on which SHIM is designed to operate 
as it has a strong band of Galactic cirrus. Cirrus structure like this 
is not typical of the HerMES fields and, in fact, the FLS field has 
been excluded from statistical studies for HerMES due to the cirrus 
contamination. While all of the HerMES maps are visually flat and 
artefact free, a program of verification, described in the sections be- 
low, was performed to determine on what spatial and point source 
flux scales these maps are reliable for use in scientific analysis. 



4 SYSTEMATICS CONTROL 

We describe below the various tests we apply to test for systematic 
effects in the SHIM maps. This includes a correction of the abso- 
lute astrometry (Section [4T| ), detailed analysis of the noise proper- 
ties of the maps based on jackknife (difference) maps (Section [4!2] 
& Section [43] >, measurement of the map transfer function in both 
deep and wide maps (Section [44) based on simulated beam-scale 
and large-scale maps and finally a characterization and bound on 
the effect of a small, correlated relative pointing error seen in the 
deepest maps (Section |43] ). 

4.1 Absolute Astrometry 

The absolute astrometry of the data is measured by stacking the 
SPIRE maps at the positions of MIPS 24 \im and radio sources. 
24 \im source catalogs come from the Spitzer Multi-Band Imag- 
ing Photometer (MIPS) ( [Werner et aL[ ( [2504} , |Rie'ke et al.|p004] )). 
These include the publicly available GOODS catalog and cata- 
logs fr om [Shupe et al.|p058] ) (Lockman- SWIRE) and |Fadda et~aT] 
(2004) (FLS). Comparisons were also made to radio source cata- 
logs in all fields including |Owen & Morrison (2008 ) (Lockman- 
North), though the corrections to the astrometry are based on the 
24 \im catalogs. Fig. [2] shows the results of the stacks along with 
the centroids and FWHM ellipses resulting from a 2d Gaussian 
fit to the stack where a 2d fit is chosen due to the known, small 
mean ellipticities of the SPIRE beams of (7 %, 12 %, 9 %) cfGrTF] 
|fin et al.|2010| >. The Figure also shows that while the 24 \im stacks 
are slightly broader than the nominal beam sizes, particularly at 
500 \im, the radio stacks are well described by the nominal beams 
sizes of 18.1, 25.2 and 36.6arcsec ( Swinyard et al.|2010| . 

The absolute pointing accuracy of the MIPS catalog sources 
is better than one arcsec, as determined from matching 2MASS 
sources in the SWIRE MIPS fields. Radio stacks are consistent 
(within 0.5 arcsec) with the MIPS stacks in Lockman-SWIRE and 
the other HerMES fields and an individual overall correction to the 
absolute astrometry, on the order of a few arc seconds, was made to 
each SPIRE map based on the 24 |^m stacks. 

In all cases, the number of sources used in the MIPS stacks far 
exceeds the number of sources for the radio (by factors of 10, 4, 30, 
8, and 100, for the five fields). We correct the vl.O maps based on 
the higher-SNR MIPS stacks (particularly for Abell 2218, where 
there are only 10 radio sources in the catalog). The offsets are not 
consistent from field-to-field since the pointing offset is known to 
be fixed in telescope coordinates and not celestial coordinates, so 
depending on the position angle of a given field, the offsets will 
be different in celestial coordinates. These offsets have now been 
corrected in the hipe pipeline and should no longer be necessary at 
the map-making stage. 

4.2 Jackknife Tests 

We divide the data into two halves in various ways to examine the 
noise properties. Figure [3] shows the distribution of signal-to-noise 
ratio from the three types of jackknife difference maps, 'orienta- 
tion', 'time' and 'detector', for each band. 'Orientation' splits the 
data by the position angle of the scans across the field into (nearly) 
orthogonal sets. 'Time' splits the data into halves based on obser- 
vation time. 'Detector' divides the focal plane at each band into two 
sets of detectors, each of which approximately covers the full spa- 
tial extent of the focal plane. The sub-maps are made with the map- 
ping parameters {p n ds , g d and w ds ) determined from the full data set. 
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Figure 2. The GOODS-N SHIM maps stacked at the location of 24 |Lim 
(upper panels) and radio sources (lower panels). Signal-to-noise ratio im- 
ages and radial profiles at 250, 350 and 500 |j.m, from right to left. The 
FWHM ellipses corresponding to the 2-d Gaussian fits are overlaid on the 
images. The major and minor radial profiles (red dashed lines) and nominal 
beam profiles (blue solid lines) are also shown. Pixel sizes are 6, 8.33 and 
12arcsec at 250, 350 and 500 |j.m, respectively. Peak positions are used to 
determine pointing offsets and the good agreement with the nominal beam 
shapes indicates the result is not heavily biased due to clustered, confused 
sources in the stack. 
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Figure 3. Distribution of map signal-to-noise ratio measured from jackknife 
maps in the GOODS-N field. Three different cuts of the data are shown: (i) 
'orientation' compares scans across the field at one position angle compared 
to the orthogonal direction; (ii) 'time' compares the first and second halves 
of the data, split by observing time; and (iii) 'detector' compares data col- 
lected from one half of the detectors in each band (alternating detectors 
are chosen) compared to the other half. The best-fit Gaussians are shown 
as dashed red curves and the expected cr = 1 Gaussians are shown as dot- 
ted lines. Below each histogram is a plot of the histogram residuals with 
the best-fit residuals in red and the residuals from the <x = 1 histograms in 
black. 



For each band and type of jackknife, the sub-maps are subtracted 
and the distribution of pixel SNRs are examined. We see from the 
overplotted best-fit Gaussians and the residuals below each his- 
togram that they are very well described by Gaussian curves, as 
one would hope if the map noise is uncorrelated. The 1 cr Gaus- 
sian widths of these histograms are: Orientation (0.95, 0.95, 0.94), 
Time (0.91, 0.91, 0.90), Detector (0.95, 0.95, 0.89) at 250, 350 and 
500 |^m respectively. Assuming that the jackknife difference maps 
give an accurate measurement of the map noise, the slight deviation 
from the expected zero-mean, cr = 1 Gaussians (particularly notice- 
able in the the black residual curves) indicates that the noise maps 
are biased high by 5 to 10 percent in all cases. For this reason, we 
recommend using the jackknife difference maps to accurately as- 
sess the noise in the maps. 



4.3 Noise Power Spectra 

Auto power spectra of both the 'time' and 'orientation' jackknife 
difference maps were produced from the GOODS-N maps after ac- 
counting for partial sky coverage. The power spectra are shown 
in Figure [4] These power spectra were produced for the map built 
with a simple, 'naive', projection algorithm as well as with the it- 



erative SHIM algorithm. In the naive maps, the 'orientation' jack- 
knife contains more residual noise than the 'time' jackknife since 
the two different scan orientations in the 'time' map allow better 
removal of the 1 // noise. In both cases, however, the iterative map- 
maker allows us to mitigate this effect on scales larger than the 
instrument beam. The 'time' jackknife spectra are almost identical 
whether they are computed with the 'naive' or iterative map-maker 
due to the large number of scans and the two different scan ori- 
entations. For comparison, the cross- spectrum of the astrophysical 
signal in the FLS field is also shown (in black) in Figure [4] This 
spectrum is obtained by correlating two maps of the FLS field us- 
ing two independent halves of the FLS data, split in time. Since 
the noise between the two halves of the data is uncorrelated, the 
cross-spectrum shows the astrophysical power with the instrument 
noise removed. The FLS field contains more Galactic cirrus signal 
than the GOODS-N field, but the extragalactic background, domi- 
nant at k larger than 0. 1 arcmin -1 , should be statistically equivalent. 
The FLS cross-spectrum is therefore a good approximation of the 
astrophysical power in GOODS-N and Fig. [4] shows that the as- 
trophysical power is dominated by the instrument noise on scales 
smaller than the beam FWHM in each band." 

The instrumental noise power spectrum amplitude, compared 
in the same region, is 10 to 100 times lower than the astrophysical 
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Figure 4. Auto power spectra of the difference jackknife maps of the 
GOODS-N field, separating the scans both by time (blue lines) and by scan 
orientation (red lines), for each band (250 to 500 \xm from top to bottom) 
with both naive (dashed lines) and SHIM maps (solid lines). These show 
that even in maps with only a single scan orientation, SHIM removes low 
frequency noise on scales below 1 degree. For comparison, the solid black 
lines represent the cross-spectrum of the astrophysical signal as measured 
in the FLS field. 



signal on scales larger than the beam. On scales larger than 10 ar- 
cminutes, there does appear to be some additional power left com- 
pared to a pure white noise spectrum. This extra power may origi- 
nate from residual temperature drifts which are currently only mod- 
eled by 2nd order polynomials. In order to trace and bound the ef- 
fect of this excess power on specific scales, in the next sub-section, 
we describe the use of simulations to determine point source fi- 
delity and to determine the SHIM transfer function (Section [4. 4\ . 
We investigate a small but correlated relative pointing error (Sec- 
tion |4.5| > and we also comment in Section [5] on low level residual 
effects due to timeline processing of the SPIRE data that may re- 
main at this early stage of the Herschel mission. 



4.4 Map Transfer Function 

The map transfer function for the SHIM SDP maps, i.e. the correc- 
tion in Fourier space for the effects of the mapper based on the ratio 
of the Fourier transforms of the simulated maps and the input skies 
for that simulation, has been measured using simulated datasets in- 
cluding beam- scale-only power, and separately, including Gaussian 
fluctuations on all scales. A deviation of the transfer function from 
unity indicates power either added or suppressed in the maps by 
the mapper itself. Where deviations from unity are small (< + 10- 
15 %) measurements of angular power in a map must be corrected 
with this transfer function to recover the true power on the sky. 




10 15 

[mJy] 



Figure 5. One of the 50 simulated noiseless realizations of the BLAST sky 
are shown in both input (left) and SHIM output (right) at 250 \xm. The input 
map is simply a simulated point source field convolved with a Gaussian 
beam of FHWM equal to that of the SPIRE beam at 250 fim. The output 
map shows the characteristic scan pattern of the SPIRE bolometer arrays as 
the input map has been resampled into timelines with the exact scan pattern 
of the real observed map. The resampled timelines are then run through the 
SHIM to produce the map on the right. 



Larger deviations indicate that the maps are no longer reliable at 
the corresponding angular scale. 

Using the best fit galaxy counts model from BLAST (the Bal- 
loon Borne Large Area Submillimeter Telescope, Patanchon et al.| 
2009 ), fifty noiseless realizations of the 250, 350 and 500 \im sky 
were created. While improved counts from |Qliver et al.| ([2010a) 
at bright fluxes and Glenn et al. (2010) at faint fluxes now super- 
sede the BLAST counts, these counts provide a reasonable model 
of the SPIRE sky and were readily available for SPIRE SDP work. 
Simulated skies were resampled into timelines and fed, in the same 
way the real data, into the SHIM iterative mapping pipeline. In- 
put and output maps for one of the 50 realizations of the simulated 
sky are shown in Figure [5] Figure [6] shows the histogram of the 
average of 50 difference maps made by subtracting the simulated 
input point source sky (left panel of Figure [5} and the SHIM map 
(right panel of Figure [5}. A Gaussian with a width corresponding 
to the (lcr) instrument noise in the deep GOODS-N map (<T inst = 
1.7 mJy beam -1 ) is shown overplotted in red. Also shown as a dot- 
ted line is a Gaussian with width corre sponding to the (lcr ) con- 
fusion noise (<r con f = 5.8 mJy beam -1 , |Nguyen et al. < 2010 )). The 
interquartile range of the histogram is 0.4 mJy beam. Since this 
includes any errors in rendering the point sources and any other 
signal introduced by the mapper, this range gives an indication of 
an upper limit on the spread in point source errors. For a source at 
the SPIRE 250 \im confusion limit (3cr con f « 18mJy), this corre- 
sponds to a 2.2 % error. The current SPIRE flux calibration accu- 
racy is quoted as +15 % ( Swinyard et al.|20T0| ). Therefore, errors 
due to the map-maker at the 2 % level are negligible in relation to 
the current calibration uncertainties. 

Testing for mapping effects on large-scale structures has been 
performed using 100 simulated realizations of beam-convolved 
Gaussian sky fluctuations, as presented in Figure [7] Shown are the 
auto- spectrum transfer functions at all three SPIRE wavelengths 
for two different fields, GOODS-N and Lockman- SWIRE, with re- 
spective map widths of 30 arcmin and 4 deg and scan speeds of 
30 and 60 arcsec/s. As in the above analysis, simulated maps are 
resampled into timelines which are then fed to SHIM. To inter- 
pret these transfer functions, the dot-dashed green lines show the 
expected transfer function due to the pixel size effect, described 
by the function (sinc(& x # pix /2) x smc(& y #p ix /2)) 2 , a necessary con- 
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Figure 6. Pixel histogram of the average of the noiseless, simulated input 
maps minus the resulting SHIM output maps. Overplotted in red is a Gaus- 
sian with cr = 1.7 mJy beam -1 , the width set by the instrument noise in 
the deep GOODS-N map. Also shown (dotted) is a Gaussian showing the 
confusion noise with cr - 5.8 mJy beam -1 . The interquartile range of the 
histogram is 0.4 mJy beam -1 . 



sequence of using a finite pixel size. This effect is detectable in 
the transfer function since the input simulations have a much finer 
(2 arcsec) resolution. Apart from this pixelization effect, which 
must be corrected for in measurements of the angular power in the 
SPIRE maps, the mapper transfer function is approximately unity 
from scales from around quarter of the size of the survey (roughly 
0.02 and 0.1 arcmin -1 for Lockman-SWIRE and GOODS-N re- 
spectively) to the scale of the FWHM for the relevant band (1.5 to 
3 arcmin -1 ). 

At angular scales just above the beam FWHM in each band, 
this transfer function fails. While primary beam effects largely can- 
cel in this ratio, since both the input simulated maps and the output 
SHIM maps have the same beam, the effect of the mapper slightly 
alters the beam widths. As shown above, this only amounts to a 
small error in the total point source fluxes, but a tiny amount of ex- 
cess power at these scales, where the input simulations have zero 
power, creates an explosion of the transfer function. As demon- 
strated in Amblard et al. (2010), this transfer function correction is 
no longer viable for k > 1, 2 & 3 arcmin -1 at 250, 350 and 500 \im, 
respectively. |Glenn et aL] ( |2010| ) detail the separate set of simula- 
tions required for their study of the beam scale power in the Her- 
MES maps. The iterative mapper seems to create 20 % extra power 
at about half the size of the map (here 0.01 and 0.05 arcmin -1 ). 
This effect is dependent on the scan length, and could be an effect 
of the choice of polynomial order. Current measurements on these 
scales should correct for this effect and this mapping artefact will 
be addressed in future data releases. SHIM is, however, very clean 
at most angular scales under the realistic conditions for which we 
tested it and we note that the map transfer function does not appear 
to depend on the scan speed. 



4.5 Relative Pointing 

In the deepest maps, we see a small but correlated pointing error. 
This is a correlated shift of the point source positions, generally on 
the order of 1 arcsec, with maximum errors as large as 3 arcsec, 
the detectability of which depends on the depth of the observation. 
Figure [8] shows a portion of the GOODS-N jackknife maps. The 
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Figure 7. Transfer function, T(k), of the simulated large-scale structures in 
the wide Lockman-SWIRE (top) map, observed at 60 arcsec/s and the deep 
GOODS-N (bottom) map, observed at the nominal scan speed of 30 arc- 
sec/s. The transfer functions are well described by the pixel window func- 
tions (green dot-dashed lines) from 0.02 arcmin -1 to about 1 to 3 arcmin -1 . 
The transfer functions show about 20 per cent extra power on the scale of 
half size of the map (0.01 arcmin -1 (top) and 0.07 arcmin -1 (bottom)). 



left panel is the 'time' jackknife in which the residual noise ap- 
pears white. The right panel shows the 'orientation' jackknife in 
which we see small, correlated residuals rising above the noise. As 
seen in the nearly Gaussian histograms in Figure [2] the effect on 
the overall map statistics is small. The artefact has not been fully 
characterized at this early stage of the Herschel mission, however, 
and we urge caution for those performing detailed statistical analy- 
ses of the SPIRE maps, whether they are the maps described in this 
work or maps produced independently by HIPE or another, custom 
map-maker. 

In order to determine the angular scales on which this issue 
is important and to bound a worst-case limit, we developed a sim- 
ple field 'warping' model which creates correlated pointing errors 
throughout a simulated map. The warping effect is visually rem- 
iniscent of a much more pronounced warping seen in early data 
due to the timing offset between the SPIRE and spacecraft clocks 
discussed in Section [2] This issue of clock drift is now well un- 
derstood and a permanent solution has been implemented in the 
HIPE pipeline. However, while the underlying cause of the remain- 
ing warp is likely due to another effect and is so far not well un- 
derstood, a model for the field warping that produces an error with 
similar features to the observed error can be created easily by re- 
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Figure 8. Time (left) and Orientation (right) jackknife maps of the GOODS- 
N field smoothed with an 18 arcsec FWHM Gaussian kernal. The right 
panel shows a small residual correlated pointing error enhanced enhanced, 
for visual purposes, by the smoothing function. 

introducing a continuous clock drift into a simulated observation 
between the detector positions and the associated pointing. The 
modeled drift was chosen to create an amplitude of the errors five 
times larger than the largest shifts seen in the sky maps. Figure [9] 
shows the ratio of the power in the 250 \xm simulated map with 
warping to the power in a map from the same simulation without 
warping. The excess power seen near 4 arcmin -1 indicates that this 
effect will be seen most strongly in the map power on angular scales 
just below the beam diameter. In real data, the effect in the science 
maps is small, as seen in the jackknife power spectra in Fig. [4] 
While the white noise level is slightly higher in the 'orientation' 
jackknife maps, there is no significant bump near 4 arcmin -1 as 
one might expect given the simulated P(k) in Figure [9] This indi- 
cates that this is a minor issue, but it is one that in not yet corrected 
for in either HIPE or SHIM and should be noted. The effect of 
this shift has been simulated and bounded by those HerMES SDP 
papers that are sensitive to these kinds of issues, e.g |Glenn et"aT] 
( |2(H0t and |Amblard eTaTl ( [20T0l ). 



5 CONCLUSIONS 

The HerMES team has designed and implemented a map-making 
pipeline, called SHIM, for the production of maps from hipe- 
processed, SPIRE data observed in SPIRE Large Map or SPIRE- 
PACS Parallel modes. SHIM reproduces both point source and dif- 
fuse extended emission with the highest possible fidelity for the 
most sensitive statistical analyses in maps with pixel sizes of 6.0, 
8.3 and 12.0 arcsec in the SPIRE bands centered at 250, 350 and 
500 \im. A rigorous program of investigation into the noise prop- 
erties of the SHIM maps, their relative and absolute astrometry, the 
map transfer function and low level residual effects, shows that the 
HerMES SDP maps are reliable on scales from the beam diameter 
up to at least 1 degree. Transfer function corrections to both point 
sources and large-scale structure are at the level of < 5%. Since 
these effects are smaller than the current overall flux calibration er- 
ror of + 15 %, we have determined that the vl .0 SHIM maps are re- 
liable for large-scale structure studies of the FIR/submm Universe 
using the bulk statistical properties of the SPIRE maps on angu- 
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Figure 9. The ratio of the power in a 250 fim simulated map with a highly 
exaggerated pointing error (warping) model to the power in a map from 
the same simulation without pointing error. The warping model, meant to 
bound the effect of the pointing error seen in Figure [8] has pointing errors 
5 times larger than the largest errors seen in the true maps and cover the 
entire field. In reality, only some regions of the map are affected and the 
amplitude of the shift is less than 3 arcsec. The excess near 4 arcmin -1 , 
however, indicates the scale on which this pointing error may be important 
for statistical studies. 

lar scales from beam size up to about half the map area. Resid- 
ual mapping errors at the few percent level, such as the pointing 
drift artefact, and their effect on sensitive SDP science have been 
flagged here. The HerMES statistical analyses most at risk of be- 
ing affected by such residuals, Glenn et al. ( 2010 ) and Ambl ard et| 
[aT]([2010), have bounded the effects on their science results using 
simulations as described in those papers. The SHIM maps exhibit 
about 20 per cent excess power at scales larger than one degree and 
analysis requiring these scales should address the potential effect 
of this excess power on their results. Timeline processing modules 
in hipe that deal with temperature drift removal and cosmic ray re- 
moval are still under development. Currently, the temperature drift 
removed from the detector signal is based on 2nd order polynomial 
models of the thermistor signal. Residual temperature drift left in 
the detector signal, as well as potential glitches missed by the cur- 
rent glitch removal modules could conceivably be leaving artefacts 
in the maps which have been bounded here and are expected to be 
continually improved. As in the case of the pointing error, those sci- 
ence efforts aimed at pushing the limits of the SPIRE data should 
characterize and bound the relevant effects on their science results 
based on the SDP processing of the SPIRE data. As the Herschel 
mission moves from SDP into Routine Observations, we anticipate 
improvement in both the hipe timeline processing and map making. 
Increasing sophistication in the SHIM implementation will likely 
lead us toward to an optimum map-maker that, along with improved 
timeline processing will extend the reliability of the HerMES maps 
well beyond the limits described here. The software discussed here 
is scheduled for release by HerMES along with the DR1 and is now 
being coded by the NHSC for inclusion in HIPE. 
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